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ABSTRACT 

Recently, mid-infrared instruments have become available on several large ground-based telescopes, 
resulting in data sets with unprecedented spatial resolution at these long wavelengths. In this paper 
we examine 'ground-based-only' diagnostics, which can be used in the study of star-forming regions 
in starburst galaxies. By combining output from the stellar population synthesis code Starhurst 99 
with the photoionization code Mappings, we model stellar clusters and their surrounding interstellar 
medium, focusing on the evolution of emission lines in the N- and Q-band atmospheric windows 
(8 - 13 and 16.5 - 24.5 /Ltm respectively) and those in the near-infrared. We address the detailed 
sensitivity of various emission line diagnostics to stellar population age, mctallicity, nebular density, 
and ionization parameter. Using our model results, we analyze observations of two stellar clusters 
in the overlap region of the Antennae galaxies obtained with VLT Imager and Spectrometer for 
mid Infrared (VISIR). We find evidence for clumpy, high density, ionized gas. The two clusters are 
young (younger than 2.5 and 3 Myr respectively), the surrounding interstellar matter is dense (> 10"* 
cm~^) and can be characterized by a high ionization parameter (logU > -1.53). Detailed analysis 
of the mid-infrared spectral features shows that a (near-)homogeneous medium cannot account for 
the observations, and that complex structure on scales below the resolution limit, containing several 
young stellar clusters embedded in clumpy gas, is more likely. 

Subject headings: galaxies: individual (NGC4038/4039) — galaxies: starburst, star clusters — ISM: 
HII regions — infrared: ISM 
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1. INTRODUCTION 

With a star formation rate (SFR) of several tens to 
several hundreds of solar masses a year, starburst galax- 
ies form stars at a much higher rate than regular gas- 
rich galaxies. As a result, starburst galaxies empty their 
gas reservoirs quickly during a short, vigorous burst of 
star-formation. This short-lived starburst phase is gen- 
erally initiated by interactions with neighbor galaxies, 
causing distortions in the galaxy's gravitational poten- 
tial. As a consequence molecular gas piles up, and in 
the resulting concentrations of very dense gas active star 
formation takes place. Although locally rare, because of 
the limited duration of the starburst phase and the rela- 
tively low galaxy pair density in the local universe, star- 
burst galaxies are a very important piece in the puzzle of 
galaxy evolution. In a very short time the galaxy's mor- 
phology and stellar content are significantly altered, and 
the interstellar and even intergalactic medium is rapidly 
enriched. At higher redshifts starburst galaxies are a 
much more common phenomenon, and sub-mm obser- 
vations indicate that the intense outbursts of star for- 
mation in these objects might be the do minant mode of 
star formation in the high-z universe fe.g. lHeckmMlll998t 
[Blain et al. 1999; Small ct al. 2002, and many others). 
Despite numerous studies of starburst galaxies, several 
issues remain unsolved, like the exact mechanism that 
ignites and shuts off the starburst phase, and the mass 

^ Based on observations collected at the European Southern 
Observatory, Paranal, Chile, under program no. 075.B-0727(A), 
075.B-0791(A), and 69.B-0688(A) 

2 Leiden Observatory, Leiden University, PO Box 9513, 2300 RA 
Leiden, The Netherlands 

^ University of Hawaii, 2680 Woodlawn Drive, Honolulu, USA 



function of the individual stellar populations formed dur- 
ing the bursts. 

In this work we focus on the individual star-forming 
regions formed during these starbursts. With a wealth 
of young, bright stellar clusters, starburst galaxies are 
excellent objects for the study of the physical processes 
of star formation. This is why starburst galaxies play 
an important role in the debate on the universality o f 
the stellar Initial Mass Function (IMF, iKroupal [2007b . 
Since star formation takes place in dusty giant molecu- 
lar clouds, the active star- forming regions in starburst 
galaxies generally suffer from high obscuration. This 
makes the (mid-)infrared a favorable regime for the study 
of the youngest, most massive, deeply embedded stellar 
populations. Analysis of mid- infrared nebular emission 
lines of clusters in starburst galaxies have led to vari- 
ous different conclusions on the IMF of these objects. 
Contrary to what one would expect for young, massive 
stellar populations, ISO measured relatively low values 
for the [Ne III] /[Neil] ratio , indicating a soft UV radia- 
tion field. iThornlev et al.l ([2000) suggest that the soft 
radiation field is a result from the aging of stellar pop- 
ulations with a 'regular' IMF (up to Mup ~ 50 - 100 
Mq), formed in short-lived bursts of star formation (1 - 
10 Myr). Alternative explanations are a truncated IMF, 
lacking massive stars, or massive stars spending a con- 
siderable fraction of their lifetime embedded in an ul- 
tra compact (UC)H ii regi on, hidden from view even at 
mid-infrared wavelength s (iRigbv fc Riekil2004D. In re - 
lation to this discussion, iMartin-Hernandez et al.l ()2002f ) 
already pointed out the need for a more comprehensive 
study of mid- infrared diagnostics, since the neon ratio is 
not only sensitive to radiation hardness, but to mctallic- 
ity, nebular density and ionization parameter as well. 
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Over the last decades the field of mid-infrared as- 
tronomy received an enormous impulse from observa- 
tions with space telescopes; from Infrared Astronom- 
ical Satellite (IRAS) in the eighties, Infrared Space 
Observatory (ISO) in the nineties and currently from 

th e Spitzer Space Telescop e. However, as discussed 

bv lMartfn-Hernandez et alT (|2005D . the interpretation of 
large aperture mid-infrared observations obtained with 
space telescopes can be challenging, especially for ex- 
tragalactic objects. Ground-based mid-infrared astron- 
omy faces additional challenges due to the characteristics 
of the earth's atmosphere, but thanks to rapid progress 
in detector technology, several mid-infrared instruments 
have become available on large ground-based observato- 
ries in the past five years. Compared to observations 
from space, ground-based work offers the considerable 
advantage of a more detailed view of the objects of in- 
terest, due to the much higher spatial resolution that 
can be obtained. For example, the typical size of a star- 
forming region is a few to a few tens of parsecs. With its 
85 cm mirror, Spitzer can resolve details on a 50 parsec 
scale in objects out to 4 Mpc distance. With an 8 meter 
class telescope similar scales can be resolved out to 40 
Mpc, increasing the volume that can be probed at this 
resolution by a factor of a thousand. Unfortunately, the 
atmosphere is opaque to most infrared radiation, mak- 
ing it impossible to access the full infrared wavelength 
regime from the ground. This means that we are limited 
to a number of specific atmospheric windows of (reason- 
ably) good transmission; the L band between 3.0 and 4.0 
/im, the M band in the range of 4.6 to 5.2 /im, the N 
band from 8 to 13 /im and the Q band between 16.5 and 
24.5 /im. The aim of this paper is to explore the spectral 
diagnostic features available for the study of star- forming 
regions within these windows, combining near- and mid- 
infrared wavelengths (J, H, K, L, N and Q band, roughly 
covering 1-25 fj,m). We will address the behavior of 
near- and mid-infrared emission line ratios as a function 
of cluster age, gas pressure, metallicity, and ionization 
parameter. 

In Section 2 we discuss the spectral features observ- 
able with ground-based facilities and the codes used in 
our model efforts, the stellar population synthesis model, 
Starburst 99, and the photoionization code Mappings, as 
well as the model grid input parameters. The model re- 
sults are presented in Section 3, and compared with exist- 
ing models in Section 4. In Section 5 diagnostic diagrams 
are introduced, which are compared to mid-infrared data 
of H II regions from the literature in Section 6. In Section 
7 we present mid-infrared spectroscopic observations of 
clusters in the Antennae overlap regions. In Section 8 
these data are analyzed and our results are discussed. 
Finally, we summarize our findings in Section 9. 

2. MODELING THE EMISSION LINE SPECTRA 

Our work focuses on stellar populations created during 
vigorous bursts of star formation in starburst galaxies. 
These stellar populations are young, dust-enshrouded 
star clusters, still embedded in their natal clouds. 
To model these systems we combine t he stellar pop- 
ulatio n synthesis model, Starburst 99 (jLeitherer et al.l 
I1999D. with the photoionization code Mappings Illr 
(iDopita et al]|2000l 120021 : iGroves. DopiTa fc Sutherland! 
l2004f ). We constructed a grid covering a large parameter 



space, exploring a range of metallicities, ISM densities 
and ionization parameters for clusters of ages ranging 
from to 6 Myr. 

2.1. Spectral features in the wavelength range of the 
atmospheric windows 

In the near- and mid-infrared wavelength range a 
wealth of spectral features is available to probe the char- 
acteristics of the emitting source, in our case embed- 
ded stellar clusters. For these systems, the near-infrared 
bands are dominated by bright hydrogen and helium 
recombination lines. When the clusters reach the age 
at which the red super giant stars appear, these start 
dominating the near-infrared continuum and, in the K 
band, CO absorption bands appear. The spectral region 
around 10 /im (the N band) displays a complex combina- 
tion of silicate absorption, fine-structure emission lines 
originating from various species and emission by Poly- 
cyclic Aromatic Hydrocarbon (PAH) molecules. Long- 
wards of the N band, the mid-infrared is dominated by 
thermal continuum emission of dust, with fine-structure 
lines and PAH emission bands (though much less promi- 
nent than in the N band) superposed on it. Furthermore, 
in both the near- and mid-infrared strong molecular hy- 
drogen H2 lines can be found. In this paper we focus 
exclusively on the emission lines. 

With ground-based observations, information between 
the atmospheric windows is unobservable, for instance 
the [Neiii] line at 15.56 /xm between the N and Q 
band window. For this reason it is impossible to use 
[Ne III] /[Neil] ratio, which is a good diagnostic for the 
hardness of the radiation field. Longwards of the Q band 
window there are several other fine-structure lines unob- 
servable with ground-based facilities. Most of them are 
not straightforward to interpret: both the [Oiv] and 
the [Fe 11] line around 26 /im can be strongly affected by 
shocks. The [Sill] line at 34.82 /im is also very complex 
to analyse, since it is not only excited in H 11 regions, but 
in Photon Dominated Regions (PDR) and X-ray Domi- 
nated Regions (XDR) as well. The only line with a rela- 
tively simple excitation mechanism at these longer wave- 
lengths accessible to Spitzer is the [S ill] line at 33.48 /im. 
Combined with the [Siii] line at 18.71 /xm this emission 
line forms a very useful density diagnostic, which is un- 
fortunately not available in ground-based observations. 

We will address some of the emission lines only acces- 
sible with space-based telescopes, like the [Neiii] 15.56 
/im and the [Siii] 33.48 /xm lines, in order to test our 
model results with mid-infrared databases of starburst 
galaxies (obtained with ISO and Spitzer, large databases 
observed with ground-based mid-infrared facilities are 
not yet available in the literature). However, the main 
focus of this paper is the modeling of ground-based 
near- and mid-infrared observations. The lines we have 
available in these observations are the series of bright 
hydrogen and helium recombination lines in the near- 
infrared bands plus several fine-structure lines in the 
mid-infrared: [Ann] at 8.99 /im, [Siv] at 10.51 /im and 
[Neil] at 12.81 /tm in the N band and [Siii] at 18.71 
/im in the Q band. 

2.2. Stellar population synthesis model 

We use the recently released version [5.1) of Starburst 
99 to model the spectral energy distributions (SEDs) of 
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Fig. 1. — FUV spectra as computed by Starburst 99 as a function of starburst age for upper: 0.4Zq, middle: IZq and lower: 2Zq. 
Ionization thresholds for a number of relevant elements are shown. 



stellar clusters. The code constructs these SEDs by sum- 
ming the luminosity-weighted spectra of the individual 
stars according to the chosen IMF, following theoretical 
stellar evolutionary tracks for each star from Zero Age 
Main Sequence (ZAMS) to its final stage. As we only 
study the youngest phases of star formation, the ionizing 
spectrum is dominated by O and B stars. The Geneva 
evolutionary tracks, which are optimized for modeling 
high-mass stars, are used to calculate the SEDs, ap- 
ply ing the high mass l oss rat e models as recommended 
by iMaeder fc Mevne"3 (|1994f ). Kurucz's stellar atmo- 
sphere models as compile d by Lejeune (jKurucz et al.l 
Il992t iLejeune et al.l 119971 ) are applied for stars with- 
out strong winds, for which the static, plane-parallel 
atmospheres are a good approx imation. For stars with 
strong winds ISmith et all (j2002f ) implemented non-LTE, 
line -blanketed model atm ospheres for Wolf-Rayet (W- 
R. iHillier fc Milled fl99l and O stars (jPauldrach et al.l 
l200l[ ). The switch between extended and plane-parallel 
atmospheres is as described in Ceitherer fc Hcck manl 
(11995). A detailed description of th e changes of version 
5. co mpared to 4-0 can be found iu iVazquez fc Leithereil 
()2005f ). Version 5.1 is similar to 5.0, with updates in the 
density and supernova routines^. 

We created three sets of models to cover a range in 
metallicities, O.4Z0, IZq and 2Z0(Z=O.OO8, Z=0.02 and 
Z=0.04; the solar abundance set is listed in Table[T]). For 
each of these metallicities we modeled an instantaneous 
burst of star formation, in which a million solar masses 
of stars are formed following a Salpetcr IMF between 0.1 

^ This inform ation comes from the Starburst 99 website: 

|http:7/www. stsci.edu/science/starburst99/ 1 . 



and 100 Mq. To ensure that the chosen cluster mass 
of 10^ Mq does not induce an artificial upper mass cut- 
off, the spectral type distribution of the resulting stellar 
populations was examined. At Myr the highest mass 
bin (03 star) contains several stars, confirming that the 
populations fully sample all spectral types up to 100 M©. 
The spectra are calculated in 1 Myr time steps between 
and 6 Myr and additionally at 2.5 and 3.5 Myr. For 
the 2Z0 models the spectra are calculated at 7 Myr as 
well, because of the extended lifetime of the W-R phase 
at high metallicity (see Section 3.1). 

Figure [1] shows the Starburst 99 ionizing spectra (far- 
ultraviolet (FUV) between 100 and 1000 A) as a func- 
tion of age and metallicity. The hardness of the radia- 
tion decreases with age, because the highest mass stars 
evolve off the main sequence. This effect is seen most 
clearly at the shortest wavelengths, between 100 - 300 
A. Longwards of 300 A the SEDs are very similar for 
all metallicities, particularly during the first four million 
years. At the highest energies (between 100 and 300 A) 
the sub-solar metallicity spectrum is significantly harder 
than the super-solar one. This difference increases with 
age. Higher metallicity stars spend a larger fraction of 
their high energy photons ionizing the metals in their own 
atmospheres (line-blanketing). Therefore, fewer photons 
are available to ionize the surrounding matter, giving a 
softer SED. Furthermore, low metallicity stars of similar 
spectral type have a higher effective temperature than 
high metallicity stars, also contributing to a harder spec- 
trum at low metallicity. Around 3 Myr the first massive 
stars enter the W-R phase, causing temporary harden- 
ing of the SED. This manifests itself as a bump around 
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Fig. 2. — Comparison of the FUV spectra computed with the new 5.1 Starburst 99 version and the original 4-0 version. The curves show 
the v5.1 spectrum divided by the v4-0 spectrum for various ages for upper: O.4Z0, middle: IZq and lower: 2Zq. 



the [O iv] ionization threshold between 3 and 5 Myr in 
the Z© spectra (between 3 and 6 Myr for 2Zq). In the 
0.4Zo spectra the SED shows a larger number of high 
energy photons over the whole 100 - 300 A range in the 
W-R phase, making the spectra between 4 and 5 Myr 
harder than at Myr. The sharp edge at 228 A in the 
spectra at 6/7 Myr corresponds to the Heii ionization 
edge (54.42 eV). At these ages there are few photons 
available for the excitation of species with an excitation 
potential exceeding that of He ii (like [Ne ill] and [O iv] ) . 

The FUV spectra have changed significantly since the 
Starburst 99 ^.0 version. Figure [2] shows the effect of the 
implementation of sophisticated atmospheric models, in- 
cluding metal line-blanketing and non-LTE 0-star atmo- 
spheres (w^.f only treated radiative transfer for hydrogen 
and helium and assumed LTE for stars of all masses). 
The plots show the new Starburst 99 v5. 1 spectra divided 
by the output of the old version of the code {v5.1/v4-0). 
The largest differences occur shortwards of 300 A. At 
the youngest ages (0-2 Myr for 0.4Zo and Z© and at 
Myr for 2Z0) the new version produces 5 - 20% more 
Hen ionizing photons (< 300 A) than v4-0. Between 
2-6 Myr, the FUV spectra (< 300 A) of the old ver- 
sion harden where the output spectra of the new version 
show a gradual softening with age, as discussed above. 
The relative amount of high energy photons available in 
the v5.1 spectra between 2-6 Myr can be as much as 
20% (at 250 A) to 60% (at 100 A) lower compared to 
the old version. This shows that a considerable fraction 
of photons at these energies is used to ionize the stellar 
atmospheres. Around 6 or 7 Myr (depending on metal- 
licity) the old and new version give roughly identical out- 



TABLE 1 

Solar abundance & metallicity scaling 



Element Abundance'' Depletion*^ 



H 0.000 0.00 

He -1.01 0.00 

C -3.59 -0.30 

N -4.20 -0.22 

O -3.34 -0.07 

Ne -3.91 0.00 

Na -5.75 -1.00 

Mg -4.42 -0.70 

Al -5.61 -1.70 

Si -4.49 -1.00 

S -4.79 0.00 

CI -6.40 -1.00 

Ar -5.20 0.00 

Ca -5.64 -2.52 

Fe -4.55 -2.00 

Ni -5.68 -2.00 



^ All abundances are logarithmic with respect 
to Hydrogen'' Depletion given as log(X/H)gas — 
log(X/H)isM 

put spectra. To conclude, characteristics of young stel- 
lar populations and/or their surroundings derived from 
emission lines sensitive to emission in this wavelength 
range alter significantly with these recent improvements 
of the code. 

2.3. Photoionization model 

The Starburst 99 output spectra are fed into the 
photoionization code Mappings to model the combined 
SED of a stellar cluster plus the irradiated surround- 
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Fig. 3. — Example of Mappings output spectra. Solar metallicity model with a density of 10 cm ^ and an ionization parameter g of 8 



10' 



The age increases from lowest: Myr to upper: 6 Myr, in steps of 1 Myr (color version available in electronic edition). Note 



the nebular bound-free edges, the strong recombination and fine-structure emission lines and the far-infrared peak originating from heated 
dust. 



ing gaseous medium. The code combines photoioniza- 
tion theory with radiation pressure and the effects of 
dust. The cloud's ionization structure is determined by 
the ionization parameter q at the inner cloud boundary: 



q = 



Qhyc 



(1) 



QLyc being the hydrogen ionizing photon flux, R the 
radius of the inner cloud boundary, and nion the ion den- 
sity (which equals ss 1.1 times the hydrogen density nn 
if we account for helium ionization). The ionization pa- 
rameter q measures the density of ionizing photons rel- 
ative to the atomic density. It relates to the commonly 
used dimensionless ionization parameter U through U 
= 1.1 • q/c. We assume that the object's radius R is 
much larger than the thickness of the emitting shell of 
gas. In this case a plane-parallel geometry is a good ap- 
proximation, which we applied in our models. As the 
code steps through the interstellar cloud. Mappings cal- 
culates the energy balance for each thin slab by coupling 
detailed modeling of the atomic gas processes with the 
dust physics. Photoionization, coUisional ionization and 
recombination processes are taken into account, as well 
as photoelectric heating of the gas by the dust, radia- 
tion pressure acting on the dust and charge transfer re- 
actions. In this way the (FUV) photons from the input 
spectrum are reprocessed by the ISM surrounding the 
star cluster. For each given combination of input SED, 
metallicity, density and ionization parameter the pho- 
toionization code gives a unique emission line spectrum 
(see Fig. [3] and a detailed description in Section 3). 



The dust is modeled by a combination of two types 
of grains, amorphous silicates and amorphous carbon 
grains. Together these adequately reproduce both con- 
tinuum emis sion in the infra red and extinctio n at shorter 
wavelengths (jPrain e 2003; D opita et al.ll2005f) . The gram 
sizes follow a power law distribution with a slope consis- 
tent with the grain shattering model ijJones "etal][l9961 
a = -3.3, modified to allow for a lower and upper size 
cu toff of 40 A and 16 00 A respectively (see Equation 19 
in lDopita et al.|[2005f ). Mappings divides the grains into 
80 size bins and computes the absorption, scattering, and 
photoelectric heating for each bin (in older versions of the 
code photoelectric heating was calculated for the dust as 
a whole). The results are then used to calculate the dust 
temperature, which determines the mid- and far-infrared 
continuum emission. While only two simple types of dust 
are considered in the calculation of dust extinction and 
emission, all elements with a high enough gas phase con- 
densation temperature (> few 100 K) are depleted. De- 
pletion factors are assumed to be ident ical to those in 
the local ISM (Sa vage &: Sembachl[T99l l. The depletion 
set is listed in Table [TJ 

Although PAH features have been observed routinely 
ever since their discovery two decades ago, exact knowl- 
edge on PAH shapes, size distributions and their inter- 
play with the environment still remains too poor for 
PAHs to be modeled adequately. In Mappings PAHs 
are represented by a single type, coronene (C24H12), for 
which the abundance is determined by the amount of 
carbon depleted onto PAH molecules (fraction here cho- 
sen to be 0.05). For this one specific PAH molecule ab- 
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sorption, emission, and photoelectric processes are cal- 
culated in detail to represent the PAH infrared spectrum 
(jPopita et al.ll2005l and references therein). 

2.4. Parameter space covered by the model grid 

The model input parameters for Mappings were chosen 
to mimic the wide spectrum of observed nebular condi- 
tions in Hii regions, covering low-density to UCHii re- 
gions. For a range in ages (0 - 6/7 Myr) and metallic- 
ities (0.4Zq, Zq and 2Zq) we computed plane-parallel, 
isobaric, dusty models with gas densities increasing from 
10^ to 10^ cm . Likewise an appropriate range in values 
for the the ionization parameter q (see Eq. [1]) was cho- 
sen. Typical values for the ionization parameter in local 
starbu rst galaxies span —3 < logU < —1.5 (Fig. 10 in 
iRigby fc Riek6,2004 who adopt a value for logU of —2.3 
in their models). Therefore the Starburst 99 spectra are 
scaled to obtain a fixed ionization parameter q at the in- 
ner cloud boundary at a value of 2 • 10^, 4 • 10^, 8 • 10^, 
1.6 • 10®, 4 • 10® and 8 • 10® cm s^^ (corresponding to 
logU « -3.13, -2.83, -2.53, -2.23, -1.83 and -1.53). 
Scaling the Starburst 99 spectra essentially means scaling 
the cluster mass. We will not address the cluster masses 
in present work, but we will use the mass information 
in future work on near-IR spectra of the same objects 
(Snijders & Van der Werf, in preparation). 

3. MODEL RESULTS 

Fig. [3] shows an example of the Mappings output spec- 
tra. The UV part of the spectrum originates from mas- 
sive stars at all times. The visible and near-infrared are 
dominated by nebular emission up to 4 Myr, showing 
clear bound-free edges. Beyond 4 Myr red super giants 
take over the near-infrared continuum. Additionally, 
prominent hydrogen and helium recombination lines are 
found in the visible and near-infrared. The mid- and far- 
infrared continuum originates from heated dust, showing 
the characteristic far-infrared peak. In the mid-infrared 
N band around 10 /im the SED is dominated by a com- 
bination of fine-structure lines, silicate absorption and 
PAH features. 

The peak of the hot dust continuum is seen to shift 
towards shorter wavelengths, indicating higher dust tem- 
peratures, with increasing age (and thus softer ionizing 
spectrum). This is caused by the way we define the Map- 
pings input parameters. The Starburst 99 input spectra, 
pressure and temperature are scaled to keep the density 
and ionization parameter at certain fixed values. A softer 
input spectrum requires a larger scaling factor to obtain 
the same amount of ionizing photons. This results auto- 
matically in a larger number of photons in the soft-UV 
and visible regime with increasing age. These photons 
do not affect the mid-infrared fine-structure lines, since 
these lines are exclusively excited by high energy UV pho- 
tons. The dust however is heated by both. Altogether, 
for a fixed value of the ionization parameter, there are 
more photons available at higher age to heat the dust, 
causing the far-infrared peak to shift to shorter wave- 
lengths. 

3.1. Age evolution of line ratios: mid-infrared 

Ratios of lines with different excitation potentials can 
be used to measure the temperature of the radiation field. 
Unfortunately the interpretation of these line ratios is 



complicated by the fact that they are sensitive to metal- 
licity of the ionizing cluster, ISM density and ionization 
parameter as well. This effect is demonstrated in Fig- 
ures |4] and [5l in which the evolution of various emission 
line ratios with cluster age is shown. Line ratios involv- 
ing emission lines that have (almost) identical excitation 
potentials, but different critical densities, are useful as 
density probes. 

To be able to compare our model results with 
models from the literature we discuss the temperature- 
sensitive ratios [Neiii]i5.56pm/[Neii]i2.8iMm and 
[S iv]io.5i;^m/[S Iii]i8.7i/im and the density-sensitive ratio 
[Siii]33.48Mni/ [S Iii]i8.7i/im, even though some of these 
ratios involve lines observable from space only. If we 
focus on the emission lines that are observable from 
the ground, we are limited to [AriiiJg.ggMm, [Sivjio.si^m 
and [Neii]i2.8i/im in the N band and [S Iii]i8.7i^m 
in the Q band. So, in addition we examine the 
'ground-based' line ratios temperature-sensitive ratio 
[Siv]io.5i;im/[Ariii]8.99^tm and the density-sensitive 

ratio [SlIl]i8.71^«n/[NeIl]i2.81^«n- 

The line ratios are later combined to form diagnostic 
ratio-ratio plots (see Section 5). First, [Siv]/[Siii] is 
plotted versus [Ne iii]/[Ne ii] to construct the commonly 
used diagnostic diagram sensitive to the hardness of the 
ionizing radiation. Secondly, the emission lines that 
are accessible from the ground are combined to form 
an extinction independent diagnostic diagram, plotting 
[Siv1/[Ariii] versus [S lll]/[Neii]. Adopting a standard 
iDraind ()1989t ) extinction curve, the intrinsic values of 
[Siv]/[Arili] and [S iii]/[Ncii] differ from the observed 
ratios by less than 3%, which is much less than the cal- 
ibration uncertainties of our observations. These ratios 
are thus essentially unaffected by reddening. 

Note that all mid-infrared fine-structure lines under 
study here are generated by a-elements (neon, argon 
and sulphur), so abundance ratios are expected to be 
approximately fixed. All differences in the line ratio 
behavior with metallicity originate from other processes. 

As the stellar population ages and massive stars evolve 
off the main sequence, the FUV spectrum softens causing 
ratios of high to low ionization lines like [Ne iii]/[Ne ii] to 
drop. The most massive star present at 2.5 Myr is an 
05.5 star. At later ages the line ratios show an upturn 
corresponding to the appearance of W-R stars, which 
dominate the radiation field until the age of 6 Myr (see 
Figs. H and E]). 

High metallicity stars have softer FUV spectra than 
low metallicity stars of similar mass due to line- 
blanketing within the stellar atmospheres. This explains 
the offset towards lower values in the [Ne iii]/[Ne ii] ratio 
with increasing metallicity as well as the steeper drop in 
the ratio in the first 2-3 Myr (see upper panel in Fig.|4]). 
After exhaustion of hydrogen burning in the core, stellar 
atmosphere temperatures can drop to low enough val- 
ues to initiate the W-R phase of high mass loss. Apart 
from having softer spectra initially, metal-rich stars can 
cool more efficiently, causing their W-R phase to start 
earlier than that of metal-poor stars. The first W-R 
stars appear at 2 Myr for 2Zo, 2.5 Myr for ZQand 3 
Myr for 0.4Zo. Furthermore, the phase in which W- 
R stars dominate the FUV spectrum lasts longer in the 
high metallicity case. Since metals are more abundant 



Mid-IR diagnostics 



7 



in their stellar atmospheres, lower mass metal-rich stars 
can generate enough radiation pressure to enter the W- 
R phase, where metal-poor stars of similar mass cannot. 
This effect extends the W-R phase to higher ages in the 
2Zq population, which is apparent from high values for 
[Ne III] /[Neil] up to 6 Myr. 

Given the values of the critical densities (1.8 • 10^ 
and 6.1 • 10*^ cm'^ for [Ne liljis.se^m and [Ne ii] 12.81^111 
respectively; all va l ues fo r critical densities are ob- 
tained from Ticlcnsl l2005t ) the [Neiii]/[Neii] is mostly 
sensitive to densities between 10'* and 10® cm^'^. The 
ratio decreases with approximately a factor of three if 
density increases from 10^ to 10® cm~^. Furthermore, 
the [Ne III] /[Neil] ratio is sensitive to the ionization 
parameter, because a larger value of q means that 
more ionizing photons are available per atom. The 
[Ne III] /[Neil] ratio increases more than an order of 
magnitude when q increases from 2 • 10^ to 8 • 10* cm s~^. 

The evolution of [Siv]/[Siii] shows similar behavior 
as [Ne III] /[Neil], with subtle differences due to (1) 
the specific sensitivity of the various individual lines 
to the spectral hardness of the radiation and (2) the 
density of the radiating medium (lower panel in Fig. 3]) . 
With critical densities of 3.7 • lO'^ cm~^ for [Siv] and 
1.0 • 10"^ cm^^ for [Siii], [Siv]/[Siii] is sensitive to 
slightly lower densities than [Neiii]/[Neii]. The density 
dependence is reversed, showing an increasing line ratio 
with increasing density. The sensitivity to q is stronger 
for [Siv] /[Siii] than for [Neiii]/[Neii]. As q changes 
from 2 • lOMo 8 • 10* cm s"* the [Siv]/[Siii] line 
ratio increases almost two orders of magnitude . 

Inspecting the evolution of [S iv]/[Ariii], one again 
recognizes the same shape as for [Neiii]/[Neii] and 
[Siv] /[Siii] (see upper panel Fig. E]). [Siv] /[Ann] is 
equally sensitive to q as [Siv] /[Sin] and somewhat 
more sensitive to density, because the difference 
in critical densities between [Siv] and [Arm] is 
larger than between [Siv] and [Siii] (3.7 • 10** and 
3.7 • 10^ cm~^ for [Siv] and [Ann] respectively). 

Proceeding to the density-sensitive line ratios, we note 
that [S iii]33.48^m/[S Ili]i8.7i;xm is most commonly used 
in the mid-infrared as a density estimator. Having the 
same excitation potential, the ratio of the double-ionized 
sulphur lines is mainly affected by the difference in col- 
lisional de-excitation rate. With critical densities of 1.2 
• 10^ and 1.0 • 10** cm^^ ^j. ^-^^^ jgjjj] 33 43 ^.j^g 

18.71 /im line respectively, the [Siii]/[Siii] ratio is most 
sensitive in the density range between 10^ and 10** cm~^, 
as can be seen in Fig. [51 Note that the sulphur line ra- 
tio changes lightly with q. This effect is caused by the 
fact that critical densities are mildly sensitive to electron 
temperature. So, in general, density-sensitive line ratios 
are weakly dependent on electron temperature as well; 
lines with a higher critical density favor a higher tem- 
perature. Nebulae characterized by high ionization pa- 
rameters have high average electron temperatures. This 
effect causes the decrease of [Sin] /[Sin] towards lower 
values with increasing ionization parameter. 

The excitation potentials of [Sin] and [Neii] (23.34 eV 
and 21.56 eV respectively) are close enough to make the 
ratio almost insensitive to the hardness of the radiation 



field. The substantially different critical densities, 6.1 • 
10^ cm^^ for [Nell] and 1.0 • lO"' cm"^ for [Slll], make 
[SiiiJ/[Neii] very sensitive to the ISM density in the 10'' 
- 10 cm~^ range (lower panel in Fig. [5]). Being sensitive 
to a different density regime [Siii]/[Neii] is not so much 
a replacement for [Siii]/[Siii] as a density probe, but 
nicely complementary. 

3.2. Age evolution of line ratios: near-infrared 

If we expand towards shorter wavelengths and include 
the near-infrared regime in our analysis, several other 
useful diagnostics become available. For one, the near- 
infrared is rich in hydrogen and helium recombination 
lines. As with the mid-infrared fine-structure lines these 
near-infrared lines can be used as a thermometer. In a 
soft radiation field, there are many more photons avail- 
able for the excitation of hydrogen, which has the lowest 
potential (13.6 eV for hydrogen versus 24.6 eV for he- 
lium) , resulting in weak He i recombination lines relative 
to the Hi lines. With an increasingly hard radiation field 
the helium Stromgren sphere gradually fills that of hy- 
drogen, observable by a higher ratio, until it finally sat- 
ur ates when the two St romgren spheres coincide (Fig. 5 
in lThornlev et al.ll2000( l. 

He i2.06Aim/Br7 has been shown to be an unreliable 
diagnos tic, or at least one that is very complicated to in- 
terpret (jligbxARiiHIlQOl E^msdene^ This 
is because the Hei2.o6Aim 2*P - 2*S line does not purely 
arise from a recombination cascade. The population of 
the 2*P level can be affected both by coUisional pump- 
ing from the metastable 2*S level as well as resonance 
scattering of the Hei 2*P - 2*S transition at 584 A. This 
effect makes the helium line very sensitive to nebular 
conditions. The He li.To^m originates almost exclusively 
from recombination cascade, making a Hei/Hi ratio in- 
volving this line a much cleaner diagnostic. Heii.7o^m 
and H i(10~4) (Brackett 10, BrlO) are close in wavelength 
and Heii.7o^m/BrlO is therefore hardly affected by red- 
dening. This ratio has been extensively tested as stellar 
temperature diagnostic and found to show good agree- 
ment with other temperature indicators (Rigbv & Ricki 
120041) . Unfortunately, the He ii.70/jm line is relatively 
faint, less than one-tenth of the Br7 intensity. This 
means it can usually only be detected in nearby objects 
with very bright recombination lines. 

Mappings does not predict the intensity of ei- 
ther the He I 2.06 or the 1.70 /um line. To obtain 
the He I 1.70 /im line strength we therefore scale 
the He I 4471 A line, which originates from the 
same upper level. Under the assumption of case 
B recombination and for the full range of tempera- 
tures and densities examined here the 4471 A lin e 
has to be scaled by 6.6 • 10~^ ([Beniamin et al.iri999f) . 

The age evolution tracks of Heii.7o^im/BrlO are 
flat at the youngest ages (for O.4Z0 and Zq), where 
the ratio is saturated and the Stromgren spheres are 
both maximally filled (Fig. [7|). With the incoming 
radiation field softening as the cluster ages, we see the 
Heii.70;im/Brl0 ratio drop again in a similar manner 
to the mid-infrared line ratios. The age evolution 
curves altogether have a very similar shape as the 
mid-infrared curves, showing the W-R phase upturn 
and the rapid decrease when the W-R stars disappear 
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Fig. 4. — Evolution of the emission line ratios with cluster age of upper: [Nein]i5.56;im/[Nell]i2.8inm and lower: 
[S lv]io.5ipm/[S Ill]i8.7ifim- The left panels show the curves for Z = 0.4Zq, the middle panels Z = Zq and the right panels Z = 2Zq. 
The upper row shows the case of low ionization parameter, q = 2 ■ 10^ cm s~^, the lower row that of high ionization parameter, 5 = 8- 
10* cm s~^. Curves of different line style (different color in electronic edition) represent results of models with different densities, ranging 
from 10"^ cm^'^ to lO'' cm~^. Line ratio values can be found tabulated at www.ifa.hawaii.edu/''kewley/Mappings/IRdiagnostics. 
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Fig. 5.— Like Fig. m for upper: [S iv] io.5i^m/[Ar ni]8. 99^01 and lower. [S in]i8.7iMm/[No n]i2.8i^, 
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Fig. 6. — Evolution of the [S Ill]33.4g^m/[S ni]i8.7i^ni ratio with increasing density. The sohd hnes result from the models with an 
ionization parameter of 2 • 10^, 4 • lO'^, 8 • lO'^ and 1.6 • 10^ cm (top to bottom solid line, the 2 • 10^ and 4 • 10^ are indistinguishable). 
For a higher ionization parameter the curve drops towards somewhat lower [Sin]/[Sni] values (short-dashed line: 4 ■ 10® cm s~^; long- 
dashed line: 8 ■ 10** cm s~^ ; the shaded area connects the q = 8 ■ 10* cm s~^ curves of different ages, the upper curve represents the 
average values for — 4 Myr and the lower curve shows the 5 Myr values). The differences between the curves is caused by a weak 
dependence of the [Sin] /[Sin] ratio on electron temperature (see Section 3.1 for explanation). Line ratio values can be found tabulated at 
www.ifa.hawaii.edu/~kewley/Mappings/IRdiagnostics. 
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Fig. 7.— Like Fig.|4]for Heii.roMm/BrlO. 
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Fig. 8. — Age evolution of the equivalent width of the hydrogen recombination line Br7. The left panels show the curves for Z = 
0.4Zq, the middle panels Z = ZQand the right panels Z = 2Z0. The upper row shows the case for low ionization parameter, q = 2 
■ 10^ cm s~^, the lower row that of high ionization parameter, g = 8 • 10*. Curves of different line style (different color in electronic 
edition) represent results of models with different density, ranging from 10^ cm~^ to 10^ cm~^. Line ratio values can be found tabulated 
at www.ifa.hawaii.edu/''kewley/Mappings/lRdiagnostics. 
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around 6/7 Myr. The ratio is also mildly sensitive to 
density because of the presence of dust in our models. 
When the density increases the dust competes more 
effectively with the gas for FUV photons. The age 
evolution curves do not show equally smooth behavior 
as those of the previously discussed mid-infrared ratios, 
especially at the 2Zq W-R plateau between 3.5 and 6 
Myr. This is most likely the result of numerical effects 
in the radiative transfer of helium in the Mappings code. 

Another commonly used near-infrared diagnostic is the 
equivalent width (EW) of one of the hydrogen lines as age 
estimator, usually the bright ones like Pa/3 in the J band 
or Br7 in K. The equivalent width measures the strength 
of the short-lived line emission from the young, most 
massive stars relative to the continuum emission from 
lower mass stars. Under the assumption of an instanta- 
neous burst of star formation, the EW drops with time 
(EW(Br7) is shown in Fig.[5D. The Zg and 2Z0 curves 
drop monotonically, but in the low metallicity case the 
curve shows a small jump between 3.5 and 4 Myr, cor- 
related with the W-R phase. With fewer heavy elements 
in the low metallicity W-R star atmospheres, not a lot 
of ionizing photons are absorbed in the process of line- 
blanketing and most are available for the ionization of 
the surrounding matter. These ionizing photons boost 
the hydrogen line fluxes. In the (super)solar metallicity 
case, where a lot of the ionizing photons are used to ionize 
the W-R star atmospheres, this upturn in the EW of Br7 
is not observed. The decreasing value of the EW with 
increasing density reflects the dusty nature of our mod- 
els in the same way as discussed for Heii yo^m/BrlO. For 
all but the most extreme case this results in a spread of 
less than 33% in the value of EW{Br7). Only the models 
with an ionization parameter g of 8 ■ 10* in combination 
with a density nn of 10^ cm~'^ show very distinct val- 
ues. The dust becomes the dominant absorber of Lyman 
continuum photons in these most extreme models. 

4. COMPARISON WITH EXISTING MID-INFRARED 
MODELS 



Both iRigbv fc Rieki 



(I200I 



■ , _ _ _ . . and 

iMartm-Hernandez et al.l (|2005D (hereafter RR04 and 
MH05 respectively) have performed similar studies to 
this work, but covering a more limited range of parame- 
ter space, assuming average values for parameters such 
as density and ionization parameter. In both works 
the age evolution of several line ratios was calculated 
in a comparable way to that presented here, using 
Starhurst 99 in combination with the photoionization 
code Cloudy (Fcrland 2001, RR04 used Cloudy version 
94.00, MH05 version 96.00-beta 4). Both assume stellar 
populations with a Salpeter IMF between 1 Mq and 
various upper mass-cutoffs (RR04: 30, 40, 50, 60, 75 
& 100 M© and MH05: 30, 50 & 100 Mq), exploring a 
range of sub-solar and solar metallicities. The models 
of RR04 are characterized by an ionization parameter 
of log U = -2.3 {q = 1.37 • 10* cm s~^) and a density 
of 300 cm-3. MH05 fix the density at 5 • lO"* cm'^. 
The ionization parameter for their object (NGC 5253) 
is estimated to be at least log U = -0.5 {q = 8.62 ■ 10^ 
cm s""'^). 

To accurately compare our models to those of RR04 
and MH05, we ran a set of models with identical param- 



eters (solar metallicity, q — 1.37 • 10*, and uh — 300 
cm^"^ for RR04, and solar metallicity, q — 8.62 • 10^ cm 
s~^ and Uh 5 • 10^ cm""^ for MH05). The resulting ra- 
tios of various emission lines are plotted as a function 
of age in Fig. [9l As judged from the age evolution of 
Heii.70Mm/BrlO (upper row in Fig. [9]), the temperature 
evolution of the RR04 models is very similar to ours. The 
curves for this line ratio essentially coincide. The mid- 
infrared line ratios of our models are considerably dif- 
ferent. Quantitatively the overall shape of the curves is 
similar, with an initial drop and an upturn corresponding 
to the W-R phase. However, the onset of the W-R phase 
happens earlier in our models and the relative strength of 
the upturn is different; the W-R stars do not seem to pro- 
duce an equally hard radiation field in the RR04 models 
as they do in ours. Since the discrepancy between the two 
sets of models is largest during the W-R phase, we expect 
the modeling of the W-R stellar atmospheres in Starhurst 
99 to be (partly) responsible for the differences. And al- 
though both studies use the same atmospheric libraries, 
the treatment of W-R stars seem to have changed be- 
tween v4-0 (used by RR04) and v5.1 (used in our work). 
Furthermore, there is a significant difference in the mod- 
eling of the surrounding nebula. We include dust, while 
RR04 model a purely gaseous nebula. The dust com- 
petes with the gas for ionizing photons, which will sig- 
nificantly alter the ionization structure within the neb- 
ula. There are most probably computational differences 
between Mappings and Cloudy as well, and differences 
in the assumed value for the abundances, like sulphur 
over hydrogen (S/H) (the abundance set is not given in 
RR04). 

In the lower row of panels in Fig. [9] we compare our 
model results to those of MH05. The general shape 
of the age evolution curves of MH05 shows a better 
agreement with our model results than RR04. However, 
the agreement between the values of the [Neiii]/[Neii], 
[Siv]/[Neii] and [Arm] /[Neil] line ratios produced by 
our models and those of MH05 is worse, with those of 
MH05 often over an order of magnitude larger. This 
difference most probably arises from differences in the 
geometry applied. We model the gas and dust as a thin 
slab (the thickness of the slab is much smaller than the 
distance from the radiating source to the slab's inner ra- 
dius) and MH05 adopt a spherical model with a relatively 
thick shell of gas and dust around the stellar cluster (in- 
ner radius is 0.6 pc, outer radius is 0.8 pc). Because of 
this geometrical difference, the cloud in the MH05 model 
sees an intenser UV field at the inner cloud radius for 
identical cloud parameters. This results in large differ- 
ences in predicted mid-IR emission line ratios. 

5. DIAGNOSTIC DIAGRAMS 

To construct more powerful diagnostic tools, we com- 
bine various line ratios into ratio-ratio plots. Fig- 
ure [10] shows the age evolution of [Siv]/[Siii] versus 
[Ne III] /[Neil] as a function of metallicity, density and 
ionization parameter. 

Fig. [4] already showed that [Siv]/[Siii] and 
[Ne iii]/[Ne 11] both measure the hardness of the ra- 
diation in a similar way and show very comparable 
age evolution tracks. Therefore the curves run nicely 
parallel and behave very orderly in the diagnostic 
diagram. Note however that at ages > 3 Myr the 
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Fig. 9. — Comparison of results for the age evolution of various emission line ratios of our models (solid lines) with results of RR04 
(dotted line in the top and middle row of panels) and with results of MHOS (dashed line in panels in the bottom row). The solid lines (our 
model results) and the dotted curves (RR04) in the top and middle row of panels can be characterized by Z = Zq, = 300 cm~^ and q 
= 1.37 ■ 10* cm s~^. The solid lines and the dashed lines in the lowest row of panels show our model results and those of MH05 for Z = 
Zq,(J = 8.62 ■ 10^ cm s'^ and = 5 ■ lO"* cm-^. 



metallicity curves wrap around each other, introducing 
a metahicity-spectral hardness degeneracy. This effect 
is caused by the difference in sensitivity of the ratios to 
radiation hardness. At this age the evolution curves in 
Fig. [4] change rapidly and are substantially different for 
each metallicity, causing the metallicity curves in the 
diagnostics diagram in Fig. [TO] to fold. 

The age evolution curves of [Siv]/[Ariii] and 
[S III] /[Neil] differ considerably in shape (Fig. [5]), caus- 
ing the curves in the ratio-ratio plot to look much 
more complex than the ones in the diagram discussed 
in the previous paragraph. Different behavior from 
the [Ne III] /[Neil] to [Siv]/[Siii] is expected, since 
here we plot a hardness-sensitive ratio, [Siv]/[Ariii], 
versus a ratio that is mainly sensitive to density, 
[S III] /[Neil] (Fig. [HI)- The fact that [Siii]/[Neii] is al- 
most exclusively affected by density and not so much by 
radiation hardness causes the relatively small spread in 



[S III] /[Neil] values within each individual panel. As in 
the [Ne III] /[Neil] to [Siv]/[Siii] diagram the metallicity 
curves wrap around each other at ages > 3 Myr to form 
an interwoven knot of model points that are not easily 
distinguished. 

6. COMPARISON WITH PUBLISHED DATA 

To test the model results, we compare the model out- 
put with mid-infrared data of gala ctic an d extragalac- 
tic Hii regions published by Giveon ~et al.l (|2002). This 
data set consists of a sample of 149 Hii regions; 112 
galactic Hii regions, 14 Hii regions in the Large Mag- 
ellanic Cloud (LMC), 4 in the Small Magellanic Cloud 
(SMC) and 19 in M33. From the ISO archive, the au- 
thors selected all H ii regions that are bright in one of the 
neon fine-structure lines and for which both observations 
in low-resolution full spectrum mode (Short Wavelength 
Spectrometer; SWSOl) as well as in medium-resolution 
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Fig. 10. — [Slv]/[Sin] versus [Nein]/[Nen] as a function of metallicity, density and ionization parameter. Stepping down from the top 
to the bottom panels the plots show the age evolution of the ratios in 1 Myr steps. The left column shows the lowest density case, uh = 
10^ cm""^, the middle column intermediate density, nu = 10^ cm"'' and the right column shows the very dense case, n^f = 10® cm~^. 
Within every individual panel, the blue/green curves represent the different metallicities; dark blue corresponds to Z = OAZq, light blue 
to Z = Zq and green to Z = 2Zq. The dots overplotted in the metallicity curves relate to the various values for the ionization parameter 
q, with the lowest point corresponding to q = 2 • 10^ cm and the topmost point to g = 8 ■ 10® cm s~^. All individual steps in q value 
are indicated in the bottom-right panel. The points for q = 2 ■ 10^ cm s~^ are connected by a yellow line, the points for q = 8 ■ 10^ cm 
s~^ by an orange and those for 9 = 8- 10* cm s~^ by a red li ne. Note that the m etallicity curves start wrapping around each other > 3 
Myr. The black dots are measurements of galactic Hll regions IGiveon et al.ll2002IV Data points indicating upper limits are excluded from 
this and the following plot to avoid confusion. The arrow in the bottom-left panel indicates an Ay of 50. Line ratio values can be found 
tabulated at www.ifa.hawaii.edu/~kowley/Mappings/IRdiagnostics. 
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Fig. 11. — As Fig.[T0]for [Slv]/[Arni] versus [S lll]/[Ne ll]. The arrow in the bottom-left panel indicates an Ay of 100. 
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line spectroscopy mode (SWS02) were available. The 
final spectra cover the wavelength range of 2.4 - 45 
fim with a spectral resolution r = 400 (low resolution) - 
2000 (medium-resolution) . Since these observations orig- 
inate from numerous different projects, all with their own 
distinctive research goals, the final sample is very hetero- 
geneous. It contains many regular star-forming regions 
as well as very bright and/or compact Hii regions. This 
broad range of conditions tests the ability of our model 
grid to reproduce a large variety of different H ii region 
properties. 

The data points of the Hii regions are plotted in the 
diagnostic diagrams presented in the previous section 
(Figs. [Tni and [TT]). In Fig. [TO] we see that the models 
and data agree very well in general, occupying the same 
range in the ratio-ratio values. The data points follow 
a distribution that is somewhat tilted compared to the 
lines of equal ionization parameter, most probably re- 
flecting morphological differences between the individual 
regions, resulting in variations in the ionization parame- 
ter. The horizontal spread exhibits a range in metallici- 
ties and densities. In general, the highest density models 
do not do as well explaining the data as the low and 
intermediate density models do. Furthermore, the very 
youngest age of Myr can be ruled out for most data 
points. On average the H ii data indicate an age for the 
star-forming regions > 1 Myr and an intermediate ion- 
ization parameter (4 • 10^ cm s""'^ ^ q < 1.6 • 10^ cm 

Examining Fig. [TT] shows that here again the data 
points and models overlap nicely, although in this dia- 
gram there are a handful of points displaying ratio-ratio 
values that are not covered by the model grid, mainly at 
the high ionization end. It is quite clear that a density 
of 10^ cm""^ is too high for all data points (but one). 
On average the intermediate density model curves (10^ 
cm~^) cover a significant fraction of the data points, but 
all densities ranging from 10^ to 10^ cm~^ are capable 
of reproducing the data. The mean ionization parameter 
indicated by this diagram is higher than concluded from 
the previous diagram, for most data a value of q between 
8 • 10^ and 8 • 10^ cm s"^ is required. 

7. A TEST CASE FROM THE GROUND: STELLAR 
CLUSTERS IN THE OVERLAP REGION OF THE 
ANTENNAE 

To address the diagnostics based on the spectral fea- 
tures that are accessible from the ground only, we 
use infrared data of young stellar clusters in the An- 
tennae overlap region as a test case. We combine 
near-infrared Infrared Spectrometer And Arr ay Camera 
(ISAAC) at the VLT (jMoorwood et al.l[l99l with mid- 
infrared VLT Imager and Sp ectrometer for mid Infrared 
(VISIR, iLagage et all 120041 ) spectra. 

The overlap region of the Antennae, the region of inter- 
action between the two merging spiral galaxies, is known 
from mid-infrared observations to host the most active 
sites of star formation, some severely reddened by the 
presence of massive amounts of dust. Based on ISO 
observations iMirabel et al.l (|1998[ ) concluded that 15% 
of the total mid- infrared flux at 15 /xm of the whole 
NCC 4038/4039 system originates from a deeply em- 
bedded, off-nuclear star-forming region, corresponding 
to a faint, very red cluster in the optical (source 80 in 



IWhitmore fc Schweizeij[T99^ . The clusters under study 
are this highly reddened cluster and an optically bright, 
blue cluster complex, which is the counterpart of the sec- 
ond brightest infrared sour ce in the overlap regi on. The 
source IDs are la and 2 in lSniiders et al.l ()200 6) for the 
reddened cluster and the cluster complex respectively, 
corresponding to the b right near-infrared sources IDwmc 
157 and IDwmc 136 (jBrandl et al.ll2005[) . 

7.1. Observations, data reduction and results 

The southern part of the overlap region of 
NGC4038/39 was observed with VISIR. Low-resolution 
N band spectra of the two brightest s ources in the 
mid-infrared were already published in ISnijders et al.l 
(2006). Additionally, medium-resolution spectroscopy 
around the [S III] line at 18.71 /xm was obtained for the 
same two sources. The total on-source integration time 
for these observations was 80 minutes. The same setup 
was used as with the N band observations, applying 
chopping and nodding on-slit with a 10" chopper throw. 
IID93813 was obser ved before and afte r the observation 
for flux calibration (jCohen et al.lll999l ). The data were 
reduced in the same way as the N band spectra and 
data were extracted in a 1'.'27 aperture (~ 133 pc at a 
distance of 21 Mpc^, assuming a Hubble constant of 70 
km Mpc~^; at this distance 1" corresponds to 105 
pc). 

Absolute flux calibration was obtained by normalizing 
the standard star spectrum to VISIR narrow-band fluxes. 
The accuracy of the flux calibration is 20%. Figure [T2l 
shows the resulting Q band spectra. 

For our analysis we use the N and Q band emission 
lines, and additionally the flux ratio of Heii.70;im/Brl0 
from the H band and the EW(Br7) from Kg band spectra, 
both obtained with ISAAC (Snijders & Van der Werf, in 
preparation). For the near- infrared observations the on- 
source integration time was 20 minute in each spectro- 
scopic setting. The spectra were extracted in 3'.' 5 aper- 
tures. This aperture is considerably larger than the aper- 
ture applied on the mid-infrared data. We have done so 
to capture all the hydrogen emission line flux that can 
be associated with the clusters and to compensate for a 
worse seeing in the near- compared to the mid-infrared. 

Th e data were e xtinction corrected assuming a stan- 
dard jOraini (|1989( l extinction curve. We obtained esti- 
mates for the visual extinction Ay from iMengel et al.l 
((2005). These values for the extinction were derived 
from the ratio of the near-infrared hydrogen recombi- 
nation lines Ha and Br7. For source la there was only 
one possible near- infrared counterpart, for which an Ay 
of 4.23 was determined. For source 2 we selected all 
sources within an 1'.'4 radius, resulting in four possi- 
ble near-infrared counterparts. To correct source 2 for 
extinction we use the average Ay value of these four 
sources. Ay — 0.28 (the spread in values for Ay was 
0.16 - 0.42). Table [2] lists the mid-infrared emission line 
fluxes and the EW(Br7) and He ii.7o^m/BrlO from the 
near-infrared spectra. 

The resulting mid-IR data points are compared with 
the sample of H ii regions observed in our and other 

^ Note that the distance of the Antennae is und er debate. Re- 
cently a lower distance of 13.8 ± 1.7 Mpc was found dSaviane et al.l 
120041) . which would affect the values derived here. 
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TABLE 2 

Near- and mid-infrared data of star clusters in the Antennae overlap region 



Species 


Wavelength 


Source la'' 


Source 2*^ 






measured 


ext corr"^ 


measured 


ext corr'^ 




2.166 


1.4 ± 0.2 


2.1 ± 0.2 


1.5 ± 0.2 


1.5 ± 0.2 


EW(Br7)'= 


2.166 


3.6 ± 0.4 


3.6 ± 0.4 


2.8 ± 0.3 


2.8 ± 0.3 


HeIi.70Mm/BrlO 


1.700/1.737 


0.309 ± 0.046 


0.315 ± 0.047 


0.299 ± 0.045 


0.299 ± 0.045 


[Ariii]<i 


8.99 


3.9 ± 0.9 


5.9 ± 1.4 


4.2 ± 0.7 


4.3 ± 0.7 


[S iv]<i 


10.51 


7.6 ± 0.9 


11.8 ± 1.4 


10.3 ± 2.3 


10.6 ± 2.4 


[Nell]d 


12.81 


49.6 ± 6.2 


59.8 ± 7.5 


30.2 ± 4.5 


30.6 ± 4.6 


[S iii]d 


18.71 


31.2 ± 3.8 


38.4 ± 4.7 


38.3 ± 3.9 


38.8 ± 4.0 



in IBrandT et al.l 120051) . star-forming knot B in Whit more et al. (2005)*^ Ay is 4.23 for source la 
and 0.28 for source 2"^ line fluxes in 10~^* erg s~^ cm"^*^ equivalent width in 10^^ /im 



Source la 



Source 2 




18.78 18.72 18.74 18.76 18.78 

wavelength (micron) 

Fig. 12. — Medium-resolution Q band spectroscopy of the 18.71 
/im [Sin] emission-line for left: source la, and right: source 2. 

nearby galaxies (the data set presented in Section 6). 
Fig. [T51 shows that the Antennae clusters fall within the 
general distribution of H ii regions, though right at the 
edge. 

The cluster ages are constrained using the EW(Br7) 
(see Fig. [H)) . Source la has the highest value for 
EW(Br7) and its age must be 2.5 Myr or younger. Fur- 
thermore, the metallicity of this source must be around 
solar or higher, since the OAZq cannot reproduce the ob- 
served value for EW(Br7). The age of source 2 falls in 
the range of - 3 Myr. For this source it is not possible 
to discriminate between the different metallicities. Tak- 
ing these age constraints into account, it is clear from 
the [Siv]/[Ariii] versus [S ill] /[Neil] diagnostic diagram 
(Fig. [TS]) that the highest density models, with a uh of 
10^ cm"*^, disagree with the data points. In general the 
medium density curves {tih = 10'* cm~'^) come closest 
to the data. Both clusters require a very high ionization 
parameter, the data points occupy an area in the ratio- 
ratio plot that is characterized by a value of g > 8 • 10® 
cm s^*. Further implications are discussed in Sections 
8.2 and 8.3. 



8. DISCUSSION 

8.1. Modeling infrared line emission of young 
star-forming regions 

The model grid presented in this paper covers a large 
range in metallicities, gas densities and ionization pa- 
rameters to mimic the characteristics of all kinds of star- 
forming regions. Testing the model output with galactic 
and extragalactic (UC)H ii regions shows that the models 
succeed in reproducing the mid-infrared characteristics 
of objects with a wide range of properties, from diffuse 
H II to UCH II regions as well as starburst nuclei. There 
is a very good agreement between the observed ratio- 
ratio values, like [Siv]/[Siii] versus [Neiii]/[Neii], and 
the range in values predicted by the models. 

We do not encounter some of the problems reported 
in previous work on the modeling of starburst galaxies. 
iKewlev et al.l (|200lf ) finds in the study of optical emis- 
sion lines of starbursts that the FUV-field produced by 
Starburst 99 v4-0 is too hard to be able to reproduce 
the observed line ratios (the number of photons in the 
300 - 1000 A range is too low compared to the highest 
energy photons between 100 - 300 A). As a possible solu- 
tion the implementation of line-blanketing in the atmo- 
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Fig. 13. — Comparison of the [S lv]/[Ar ill] and [S lll]/[Ne ll] ratio 
of the Antennae clusters observed with VISIR (crosses, thi s work) 
with (extra-)galactic Hll regions (dots: [Glveon et al.|[2003) . None 
of the data is extinction corrected. The Antennae data points 
include error bars. 
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Fig. 14. — This plot is identical to Fig. [8] with the data points for the Antennae clusters included. The fat black line indicates the value 
of EW(Br7) for both sources as measured from ISAAC spectra. Error bars are indicated by the shaded areas. From this plot we conclude 
that both clusters are young, source la < 2.5 Myr and source 2 < 3Myr. 



spheres of high-mass stars was suggested, which is indeed 
incorporated in the current version of Starburst 99 (v5.1). 
Nonetheless, the results of a more recent study of mid- 
infrared o bservations of starbu rst galaxies show similar 
problems (jRigbv fc Riekc 2004'). Although these authors 
use a version of Starburst 99 with line-blanketing imple- 
mented, the observed mid-infrared line ratios combined 
with He ii.70/im/BrlO can only be explained by models as- 
suming an IMF with a relatively low upper mass cutoff 
(40 - 60 Mq). In our work presented here there is no need 
to assume an unusually low M„p. As discussed in Section 
4 we expect the fact that our models include dust, while 
the Rigby & Rieke models do not, plus changes in the 
modeling of W-R stars in Starburst 99 to be responsible 
for the differences between our work and theirs. Fur- 
thermore, there may be computational differences in the 
photoionization codes used; we use Mappings and Rigby 
used Cloudy. 

In spite of the improvements, several inconsisten- 
cies remain. The average ionization parameter for 
Giveon's Hii regions derived from the [Siv]/[Siii] ver- 
sus [Neiii]/[Ncii] plot, is lower than values indicated by 
the diagram plotting [Siv]/[Ariii] versus [S ill] /[Neil] (4 
• 10^ cm s~^ < q < 1.6 ■ 10® cm s^"'^ for the first di- 
agram and 8 • 10^ cm s~^ < q < 8 ■ 10* cm s^^ for 
the latter). This discrepancy may originate from vary- 
ing relative elemental abundances and/or differences be- 
tween the detailed shape of the stellar FUV spectra in 
the models and in reality. Examining the data points in 
the [Siv]/[Ariii] versus [Siii]/[Neii] diagram (Fig. [TT]) 
suggests that the values for [Siv]/[Ariii] predicted by 
the models are too low and the strength of argon lines 



in the models is too high relative to the sulphur emis- 
sion lines. Difficulties in the modeling of sulphur hav e 
been reported befo re (jMartin-Hernandez et al.l (|2002D . 
iVerma et al.l (|2003[ )V However, in these studies on ISO 
observations of galactic H ii regions and starburst galax- 
ies the sulphur abundance was found to be low compared 
to neon and argon. 

We modeled a very simple system, a single spherical 
stellar population, formed in an instantaneous burst of 
star formation wrapped in a thin shell of gas and dust. 
We assumed a Salpeter IMF between 0.1 and 100 Mq for 
the stellar mass distribution of the clusters. Although 
the shape of the IMF is under constant debate, espe- 
cially in extreme env ironments, such as in starbursts, 
iThornlev eFall (|2000 f) showed that differ ences in metal- 
licity and ionization parameter have a much larger ef- 
fect on mid-infrared line ratios than the shape of the 
IMF. Given observational constraints, the assumption 
for the IMF upper mass cutoff is not expected to in- 
troduce a large error either. From the EW(Br7) we 
know that both sources in the Ante nnae overlap region 
are 3 Myr or younger. Furthermore. iKunze et all ()1996f ) 
derived an effective temperature of 44,000 K from the 
[Ne III] /[Neil] ratio, corresponding to the presence of a 
most massive star of 60 Mq in the present day mass 
function (PDMF). The variation of several mid- infrared 
line ratios due to a different IMF upper mass cutoff dur- 
ing the first 3 Myr is shown to be a factor of two or 
less (comparing M^p of 50 M© and 100 Mq in Fig. 3 in 
iRigbv fc Riekdl2004l ). So the choice of upper mass cutoff 
in the range constrained by the observed age and effec- 
tive temperature has a limited effect on the line ratios 
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Fig. 15. — This plot is identical to Fig. 1111 with the data points for the Antennae clusters included (crosses). The 4 Myr and 5 Myr 
panels are shaded, since we concluded from EW(Br7) that the clusters are 3 Myr or younger (source 2; source la is < 2.5 Myr). Error bars 
are given and the Antennae data points are extinction corrected as described in the text. 



compared to varying density and ionization parameter. 

The assumption of an instantaneous burst of star- 
formation does have a considerable impact on the mid- 
infrared emission hues. Because of the constant produc- 
tion of massive stars, the SED of a stellar population 
undergoing continuous star formation at a certain age is 
expected to be harder than that of a passively evolving 
stellar cluster formed in a single instantaneous burst of 
the same age. Thus under the assumption of constant 
star formation the ISM properties do not need to be as 
extreme (in density and ionization parameter) to repro- 
duce the observed line ratios. Another source of major 



uncertainty is the modeling of stellar atmospheres. Since 
the mid-infrared emission line ratios are very sensitive 
to the detailed shape of the FUV spectra, the choice of 
atmospheric models have a large effect on the analysis 
of the observations. Furthermore, we know from local 
H II regions that the morphology of star-forming regions 
is usually very complex and nothing like a simple sphere. 
Since the regions discussed here are large, with a radius 
of 40-45 pc, we expect significant clumpy structure below 
our resolution limit. This issue will be addressed in the 
next Section. Lastly we mention the presence of shocks 
as another possible excitation mechanism. Indications 
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log([SIII]/[Nell]) 

Fig. 16. — A blow-up of one individual panel of Fig. 1151 The 
metallicity curves have been extended to higher values for the ion- 
ization parameter q. The solid curves correspond to different metal- 
licities, black to Z = 0.4Zq, dark grey to Z = Zq, and light grey 
to Z = 2Zq (a color figures can be found in the electronic edition). 
The dots corresponds to g = 2 • lO'', 4 ■ lO'^, 8 ■ lO'', 1.6 ■ 10**, 4 
• 10®, 8 • 10**, 2 • 10^ and 4 • 10^ cm s'^ The grey dashed lines 
connect curves of equal q. 

for shocks being important in the Antennae overlap re- 
gion are found fro m extraordinarily bright H2 emission 
()Haas et al.ll2005l ). Part of the line emission could be 
shock excited, complicating the analysis. 

8.2. Dense star-forming regions in the Antennae 

Comparison of the model output with near- and mid- 
IR data of two clusters in the Antennae overlap region 
allow us to constrain the properties of the ISM surround- 
ing the star clusters. Fitting [Siv]/[Ariii], [Siii]/[Neii], 
EW(Br7) and Heii.7o^m/BrlO simultaneously results in 
a good fit of source 2 by a 3 Myr old stellar population 
with either solar or sub-solar metallicity, a density of Iff* 
cm~^ and an ionization parameter of 8 • 10® cm s~^. 
Given that the metallicity of most star-forming regions 
in the overlap region was found to be aro und solar from 
optical observations (jBastian et al.l [20061 ) . we prefer the 
solution with solar metallicity. 

It is evident from Fig. [15] that, unlike source 2, source 
la cannot be explained by a single component from the 
presented model grid. One possible solution to reproduce 
the observed line ratios is to extend the metallicity curves 
to higher values for the ionization parameter. Under the 
assumption that the metallicity of source la is around 
solar as well, the 2 Myr Zq curve for a density of 10^ 
cm~^ looks most promising. To estimate by how much q 
must be increased we extended this curve by adding the 
data points for a value for g of 2 • 10^ and 4 • 10^ cm s^^. 
Figure [TBI shows that the [S iv]/[Ariii] ratio saturates for 
such high ionization parameter values. The line ratios 
observed for source la can be fitted by a slightly super- 
solar metallicity, 2 Myr stellar population, embedded in 
matter with a density of 10** cm~^ that can be charac- 
terized by an ionization parameter of 4 • 10^ cm s~*. 
Compared to previo us estimates for q this is a very high 
value (see Fig. 10 in iRigbv fc Riek^ l2004l for ionization 
parameters for a collection of local starbursts. The mean 
value for logU is -2.3). However, since the star-forming 
regions in the Antennae are known to be very massive 
(~ o ne million solar masses) and bright (jMengel et al.l 
[200l . there is no reason to rule out this option before 
further analysis. A measurement of the [Nciii] line flux 
would help to obtain more robust constraints on q from 



our observations. Unfortunately it is impossible to ob- 
tain such observations with the same spatial resolution 
as our ground-based data. We note that a less extreme 
ionization parameter might be achieved by the use of 
the ult ra-compact H ii region models by 'Do pita et al. I 
()2006aD . This possibility will be explored in future work. 

Since it is known from observations that star-forming 
regions are usually made up of a clumpy, structured, 
multi-phase ISM, we consider a second option to repro- 
duce the mid- infrared characteristics of source la. We as- 
sumed a two-component medium by making linear com- 
binations of model output spectra (fi x spectrumi + f 2 x 
spectrum2, with fi and f2 varying from 0.1 to 0.9 in steps 
of 0.1 and fi + f2 = 1; we only consider combinations with 
identical metallicity). In this way one can obviously cre- 
ate a large set of possible fits to the [S iv]/[Ar ill] and 
[S III] /[Neil] values of source la. It is interesting though 
that, without exception, all these possibilities require a 
high density. At least 70% of the matter from which 
these mid-infrared lines originate must have an average 
density of > 10^ cm~^. If we fit our model results to 
the mid- and near-infrared lines simultaneously there is 
only one possible combination that can explain the data 
of source la that meets the assumption of solar metallic- 
ity: 90% of the radiating ISM can be characterized by a 
Myr stellar population with a density of 10^ cm~^ and an 
ionization parameter of 8 • 10® cm s^* and the remaining 
10% by a 1 Myr solar metallicity population with a 10^ 
cm^^ and an ionization parameter of 8 • 10® cm s~^. In 
reality a two component medium is an oversimplification 
as a model of an inhomogcneous star-forming region, but 
these results still provide a good estimate of the range of 
average properties we expect for these regions. 

Altogether combining the data and the models indi- 
cates that source la is a very young stellar population 
packed in a giant cloud of very dense gas and dust. 
During the first million years after the birth of the stel- 
lar population, one might expect the cluster still to be 
deeply buried in its natal cloud, reddened by far more 
obscuring matter than the Ay of 4.2 derived from near- 
IR data. The extinction is known to be patchy in this 
region, varying orders of magnitudes on scales compara- 
ble to our spatial resolution limit. Mid-IR narrow-band 
images revealed an ob ject close to source la suffering 
much more extinction ijSniiders et al.|[2006D . Since this 
object does not even have a counterpart detected in the 
near-IR, we concluded that the Ay must be in the order 
of 70 or larger. However, even if the extinction on source 
la is much larger than the Ay of 4.2 corrected for here, 
our results would not change, since the [S iv]/[Ar III] and 
[S III] /[No 11] diagram is relatively insensitive to redden- 
ing (see the vector corresponding to Ay = 100 in Fig.fTTjl. 

8.3. Evidence for a dense, clumpy ISM 

Comparing the average density found in this work with 
ISM densities in the Antennae overlap region from the 
literature, we find that our value is comparable to the 
density of molecula r gas derived fro m sub-mm *^C0 and 
^•^CO observations (jZhu et al.ll2003[ ). These observations 
indicate a two phase molecular gas medium, one compo- 
nent having a density n//^ of 1-8 • 10'^ cm~^ and a second 
denser component with > 3 • 10** cm~^, both average 
values for the whole overlap region. [Fe ill] emission lines 
in near-infrared spectra with comparable spatial resolu- 
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tion as our data, indicate a somewhat lower ionized gas 



density than the value we found, 10 



3.5 



10^ 



,-3 



for 



source la (under the assumptio n that uh ^ rie. So urce 
2 is not addressed in this paper, [Gilbert et all[2000l) . In 
the same work a molecular gas density of 10^ cm~'^ is 
found from the analysis of H2 line emission. 

Similar densities are locally observed only on very 
small scales in UCHll regions. Typically these star- 
forming clouds have a mol ecular gas density > 10^ cm"'^ 
within a radius < 0.5 pc (|Churc hwcll 200^. With com- 
parable average gas densities, but radii almost two orders 
of magnitude larger (the radius of source la and 2 is ap- 
proxi mately 45 pc and 40 pc respectively ISnijders et al.1 
|2006[ ). the star- forming regions in the Antennae bear 
more resemblance to extreme star-forming regions in ul- 
tra luminous infrared galaxies (ULlGs). A large fraction 
of the molecular gas in these sources is found to have 
densities larger than 10 ^ cm~^ (ISolomon et al.lll992f ). In 
a large sample of ULIGs [Solomon et al.1 ( 199^ find some- 



what lower molecular gas densities of 10 - 10 cm"'' for 
regions of 100 - 500 pc radius. In a subsequent paper, 
similar densities (2 • lO'^ - 2 • 10^ cm~^) are found for 
starburst regions in Arpl93, Mrk273 and Arp220 for re- 
gions with 68 - 150 pc radii. Gas masses for these region s 
vary from 0.6 to 1.1 • 10^ Mq ijDownes &: Solomon|[l998l ) . 

If we assume a homogeneous distribution of the gas in 
spherical clouds the gas mass for source la is ~ 1.2 • 10^ 
Mq (gas density uh of 10'* cm~'^, taking the presence of 
helium into account, and a radius of 45 pc) and ~ 8.6 • 
10^ Mq for source 2 (gas density uh of 10* cm~^ and a 
radius of 40 pc). However, the assumption of a homoge- 
neous medium is inconsistent with reality, as shown by 
the following argument. If we combine the results for the 
ionization parameter and the density derived from our 
models with the Lyman continuum emission calculated 
from the extinction corrected Br7 flux (Tabled), we can 
get an estimate for the source size using Eq. [1] Under 
the assumption of case B recombination QLyc is approx- 
imately 9.3 ± 0.9 • 10^2 g-i and 6.7 ± 0.7 • 10^^ for 
source la and 2 respectively (this is in good agreement 
with determinations of QLyc for source la by Gilbert et 
al. 2000, from near-infrared emission hydrogen lines and 
Hummel & Van der Hulst 1986, from thermal radio con- 
tinuum flux, both find 1 • lO'^^ s^*). If we focus on 
source la, and evaluate Eq.[T]with the appropriate value 
for QLyc, a density nn of 10* cm~^ and a value for g of 4 
• 10^ cm s^^, the radius of the emitting source should be 
~ 4.4 parsec (using the values for the two-component fit, 
uh = 10^ cm~'^ and q — 8 ■ 10* cm s~^, results in a radius 
of 3.1 pc). This value is an order of magnitude smaller 
than the source radius measured from mid-infrared im- 
ages, namely a half-light radius of 45 parsec. 

This discrepancy implies the need for a more realis- 
tic, physical model than a single stellar population sur- 
rounded by a single-phase gas cloud. We adopt the ge- 
ometry of a collection of N distinct H 11 regions (star- 
forming clumps) embedded in a giant molecular cloud, 
very similar to the proposed configuration in Fig. 7 of 
iForster Schreiber et al.l (j2001f ) for the nearby starburst 
M82. Observational support for this model can be found 
in mid-infrared images of source la. Although most of 
the clumpy substructure would be indistinguishable be- 
cause of the spatial resolution of our observations, the 



general idea agrees with the finding of multiple ionizing 
sources er nbedded in an envelo pe of continuum and PAH 
emission ()Sniiders et al.l [20061 ). Star clusters in cluster 
complexes ar e known t o form following a ce rtain mass 
distribution (jFalil 120061 : iDopita et all I2006U and ref- 
erences therein), so ideally one should perform Monte 
Carlo simulations of various different stellar populations 
buried in a clumpy ISM. However, for simplicity we as- 
sume N identical stellar populations. The Lyman contin- 
uum originating from each of these individual clusters is 
Qhyc/N. Combining the equation for the ionization pa- 
rameter q for each individual region with the ionization 
balance for a Stromgren sphere (Eqs. [2|and[3|), results in 
an indication of the size and number of individual H 11 re- 



gions. 



q = 



N 



N 



-7ri?|nf„„aB 



(2) 



(3) 



aB is the recombination rate coefficient, 2.6 • 10~*^ cm'^ 
s~* at Tefi — 10,000 K. R is the distance between the ra- 
diating source to the inner cloud boundary and Rg is the 
Stromgren radius. For typical H 11 regions most H 11 ra- 
diation originates from an ionized shell of gas around 
the Stromgren radius, so i? ~ i?s. Solving this set of 
equations with identical values as used above results in 
a typical radius of 1.35 parsec and a number TV of ^ 
9 clumps for source la. The Lyman continuum implies 
that within each of these clumps several hundreds of O- 
type stars are present (QLyc = 9.6 • 10^^ s~^ corresponds 
to the ionizing flux fro m ~ 225 03 stars, us ing the cal- 
ibration of O stars from [Martins et al] 12005'). The total 
ionized gas mass from all clumps added up is 2.8 • 10^ 
Mq, approximately 3.5 orders of magnitude less than un- 
der the assumption of a homogeneous cloud of gas. 

A nearby example of a system comparable to the model 
proposed here, a giant star-forming complex encompass- 
ing a handful of bright H 11 regio ns, can be found in M82 
(at 3.6 Mpc distance, see Fig.l in lSmith et"al[|2006f) . Re- 
gion A, an elongated cluster complex of ~ 50 - 100 pc, 
shows rich substructure with several tens of star-forming 
clumps. The characteristics of one of these massive stel- 
lar clusters, M82-A1, is very similar to the properties 
inferred above for the star-forming clumps: the compact 
Hii region has a radius of 4.5 ± 0.5 pc and a density 
of 1800 cm^'^. The central stellar population has an esti- 
mated age of 6.4 Myr and the total stellar mass is around 
a million solar masses. At the distance of the Antennae 
most information on the rich structure observed in M82's 
cluster complex A would be lost due to spatial resolution 
effects. It is very plausible that at a younger age the clus- 
ter complex would look very similar to one of the bright 
SSCs in the Antennae overlap region. 

Further constraints on the stellar content of the star- 
forming clumps in our model can be determined by esti- 
mating the total stellar mass for source la by comparing 
the extinction corrected Br7 flux with the Starburst 99 
output for a Salpeter IMF between 0.1 and 100 Mq. De- 
pending on age, which we constrained to be - 2.5 Myr 
for source la, the total stellar mass is between 3.2 ± 0.3 
■ 10^ Mq (0 Myr) and 4.2 ± 0.4 • 10^ Mq (2.5 Myr; at 
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the lower distance as found by ISaviane et all (|2004[ ) the 
mass would be lower by a factor of ~ 2.5). This is in good 
ag reement with the val ue of 3 million solar masses found 
bv lMengel et al. | (|2001D. but c onsiderably lower than the 
16 ■ 10^ bv [Gilbert et all jMOO). Note however that 
the studies of Gilbert and Mengel both use a lower mass 
cutoff of 1 Mq, and Gilbert adopts an age of 4 Myr. 
With a total stellar mass of 3 - 4 million solar masses, 
each individual star-forming clump would host 3 - 4 • 
10^ Mq of stars, resulting in a stellar density of 3 - 4 • 
10^ M0 pc~'^, which is equivalent to a density nn of 0.9 
- 1.2 • 10^ cm~^. This implies very high star- formation 
efficiencies if we assume these stellar populations to be 
formed from typical dense molecular clouds (cloud core 
densities nn ~ 10^ - 10 ^ cm"'^). Studies of the most 
massive galactic star-forming regions and of super star 
clusters (SSC) in various starburst galaxies found com- 
parable eq uivalent densities, as summarized in Fig. 1 by 
VTaH ([20061 ). In this Fi gure source la would end up be- 
tween the Arches cluster in the Milky Way and the cen- 
tral SSC in NGC5253. 

Although our findings for the average stellar mass den- 
sities are similar to results for other massive star-forming 
regions, there are scenarios that require less extreme 
mass densities. For a stellar mass distribution with rela- 
tively little low-mass stars compared to 'Salpeter' (e.g. a 
Miller-Scalo, Kroupa, Chabrier, or top-heavy IMF), the 
mass within the 1.35 pc radius would be less. Indications 
for the existence of an IMF dominated by massive stars 
is found in several massive star-forming regions, in the 
center of our galaxy (Nayakshin &_Sunyacv 2 005,) a,s well 
as in other starb urst galaxies ( Doane fc Mathews! Il993l : 
iRieke et al.lll993l ). Such a deviant IMF could be the re- 
sult of stellar winds and supernova explosions of the high- 
mass stars shutting off the slower formation process of 
low-mass stars by expelling the available matter from the 
star-forming region. Another option is mass segregation 
within the stellar cluster. The high-mass stars in clus- 
ters are often observed to be more centrally concentrated 
than the low-mass stars. This is possibly caused by rapid 
inward migration of mass ive stars (Freitag ct al 2006; 
iPortegies Zwart et allll999t) . The mid-infrared emission 
lines we observe, originate predominantly from the high- 
est mass stars, and the cluster core, where most of these 
massive stars reside after mass segregation, would be 
deficit in low-mass stars. This mechanism would not af- 
fect the total cluster mass, but would lower the derived 
density in the cluster core. 

Another plausible model for the morphology of a clus- 
ter complex would be several star-forming clumps em- 
bedded in a large cloud of diffuse, ionized gas. However, 
from our mid-infrared data we do not see evidence for a 
significant fraction of low density ionized gas. Gas with 
a density of 10^ cm~'^ would boost the fluxes of the emis- 
sion lines with the lowest critical densities (e.g. [S ill] at 
33.48 /im). Consequently, the presence of a large amount 
of low density ionized gas would leave its fingerprints on 
the observed line ratios, which would have been recov- 
ered as possible solutions in the two-component medium 
fits. The only possibility would be a component of low 
density gas characterized by a very high ionization pa- 
rameter (g > 5 - 10 • 10^ cm s~^). As we have seen in 
Fig. [16] the [S iv]/[Ar ill] ratio saturates at high values 
of q, causing a bend in the metallicity curves towards 



lower [S iii]/[Ne ii] values. In this way it would be pos- 
sible to have a significant contribution from low density 
gas to the observed line ratios. We consider this option 
highly unlikely though. A combination of low density and 
such high ionization parameter would not arise, since q 
is inversely proportional to and we expect the diffuse 
gas to be at a relatively high distance from the ioniz- 
ing sources on average, which would suppress q (Eq. [1]). 
Furthermore, the typical ionization parameters for star- 
forming regions in starburst galaxies derived from obser- 
vations are lower (Fig. 10 in Rigbv fc Ric k e 2004) . as well 
as the values for q predicted by models ([Dopita et al. I 
l2006aD . 

9. CONCLUSIONS 

The recent advent of the first generation of mature 
mid-infrared instruments on ground-based telescopes was 
one of the main motivations for this work. With mid- 
infrared observations being performed routinely on vari- 
ous 8-meter class telescopes nowadays, rich new data sets 
rapidly become available. A lot of work has been done 
based on space-based mid-infrared observations over the 
last decades. Because of the characteristics of the earth's 
atmosphere, ground-based work is limited to the spectral 
features observable from the atmospheric windows of rea- 
sonably good transmission. So, in this work we focus on 
the features that are available from ground-based obser- 
vations. 

We have used the stellar population synthesis model 
Starburst 99 combined with the photoionization code 
Mappings to simulate mid-infrared emission lines of star- 
forming regions in starburst galaxies. We have shown 
that the models succeed in reproducing the observed line 
ratios of H ii and starburst regions with a wide variety of 
properties, ranging from diffuse to very dense systems. 

Applying the model results in an analysis of the mid- 
infrared data of SSCs in the Antennae overlap region, 
leads to the conclusion that the average ionized gas den- 
sity in these star-forming region is very high. Source 2, 
the mid-infrared counterpart of a blue cluster complex, 
can best be modeled by a 3 Myr old stellar population 
embedded in dense matter, tih = 10^ cm~'^, that can 
be characterized by an ionization parameter of 8 • 10* 
cm s~^. The observed line ratios for source la, cor- 
responding to the highly reddened cluster WS95-80 in 
[Whitmorc & Schwcizcr (1995), can only be reproduced 
by assuming an extremely high ionization parameter q 
= 4 • 10^ cm in combination with a 2 Myr stellar 
population and 10** cni~^ density. We also considered a 
two-phase ISM by fitting linear combinations of model 
output. This resulted in a younger stellar population of 
0-1 Myr, surrounded by 10^ cm~'^ gas characterized by 
g = 8 • 10* cm s-^ 

Detailed comparison of the mid-infrared emission lines 
with the images shows us that the ISM in these two star- 
forming regions is as expected far from homogeneous. As 
a possible geometrical model we propose several individ- 
ual star-forming clumps embedded in a giant molecular 
cloud. This model is supported by the morphology of the 
source in the mid- infrared images, showing two ionizing 
sources within a large enveloping cloud of continuum 
and dust emission. The typical number of H ii regions 
within the molecular cloud of ^ 40 pc radius is 9, each 
with a radius in the order of 1.35 pc (source la). Such 
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morphology is observed in nearby starbursts as well. 
Each of these clumps would contain ~ 225 03 stars and 
have a gas density of 320 Mo/pc'^ {uh — 10^ cm~^). 
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